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The spread of the time arrows from the environment to an observed subsystem is followed within a 
harmonic model. A similarity is pointed out between irreversibility and a phase with spontaneously 
broken symmetry. The causal structure of interaction might be lost in the irreversible case, as well. 
The Closed Time Path formalism is developed for classical systems and shown to handle the time 
arrow problem in a clear and flexible manner. The quantum case is considered, as well, and the 
common origin of irreversibility and decoherence is pointed out. 



Contents 



I. Introduction 1 

II. Synopsis 2 

III. Classical harmonic systems 6 

A. System of oscillators 6 

B. Time arrow and the null-space 8 

C. Coupling of time arrows 9 

D. Continuous spectrum 10 

IV. Dynamical breakdown of reversibility and causal structure 10 

A. Spontaneous symmetry breaking 10 

B. Irreversibility 11 

C. Acausality 12 

V. Classical closed time path formalism 12 

A. Action 13 

B. System of harmonic oscillators 15 

VI. Quantum systems 16 

A. CTP formalism 17 

B. Explicit breakdown of time reversal invariance 18 

C. Unitarity and decoherence 20 

VII. Summary 21 

Acknowledgments 22 

A. CTP Green-function 22 

1. Action 23 

2. Green- function 23 

3. Finite T 24 

4. Diverging T 25 

5. Inverse Green- function 26 

References 26 



I. INTRODUCTION 



According to our daily experience the time flows "ahead" and events look differently when seen in a movie played 
backward. The time arrow which points in the direction of the flow of time seems to be well defined despite the 
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invariance of the fundamental equations of motion under time reversal, t — > —t, the weak interaction put aside. There 
are two possible levels where the orientation of time might be defined. One is General Relativity where both the 
time and its arrow are constructed. The other is in quantum mechanics or thermodynamics where quantum laws or 
statistical mechanics sets the time arrow [1-3] . The issue of the mechanical time arrow is considered in this work by 
paying special attention to the preparation of the system and the way it is observed. The information we need about 
the preparation is encoded in the boundary conditions in time, imposed on the solution of the dynamical problem 
and what is assumed about the observations is that they are made in a finite duration of time. 

The main idea, borrowed from the renormalization group method [4, 5], is that in any real experience on a dynamical 
system, called full system below, we observe a subsystem only, called briefly system from now on, and can construct 
an effective theory for this subsystem. The effective theory can in principle be obtained from the detailed dynamics of 
the full system by eliminating the non-observed degrees of freedom, called environment, by solving their equations of 
motion. One can speak of irreversibility with respect to the observed subsystem only since it refers to the dynamical 
process where the environment boundary conditions induce a time arrow for the observed subsystem of the full system 
obeying time reversal invariant dynamics. It may happen that the environment induces interactions both forward 
and backward in time and no definite system time arrow is observed, the system dynamics is acausal in such cases. 

It seems to us that the emergence of irreversibility can be understood as a dynamical symmetry breaking. This view 
is motivated by a careful separation of two scenarios, the realistic one based on observations carried out in arbitrarily 
long but finite time and another, mathematically simpler idealized situation where observations may take infinitely 
long time. Though irreversibility and acausality share the same dynamical origin their manifestations arc different. 
The former is a spontaneous symmetry breaking equipped with an order parameter and is easier to recognize. The 
observation of the latter is more involved which may explain the more abundant models and dynamical scenarios, 
developed to find the dynamical origins of irreversibility than those of acausality. 

The Closed Time Path (CTP) formalism [6] is used in this work. It has already been reinvented in different quantum 
contexts such as to derive relaxation in many-body systems [7], to develop perturbation expansion for retarded 
Green-functions [8] , to find manifestly time reversal invariant description of quantum mechanics [9] , to describe finite 
temperature effects in quantum field theory [10-12], to find the mixed state contributions to the density matrix by 
path integral [13], to follow non-equilibrium processes [14], to derive equations of motion for the expectation value 
of local operators [15, 16] and to describe scattering processes with non-equilibrium final states [17]. This scheme, 
introduced here already in classical mechanics does not give rise to new equations. However it is better suited to the 
problem considered because it separates clearly the oriented and unoriented correlations in time among the dynamical 
degrees of freedom and sources, it can deal with initial condition problem for an infinite system within the framework 
of the action principle thereby allowing the use of powerful functional methods. Furthermore, it provides a variation 
principle to derive dissipative effective equations of motion and finally it makes clear that causality is not automatic 
for infinite systems. 

We start with a succinct description of our procedure in Section II, followed by the introduction of our model in 
Section III, the presentation of irreversibility as a spontaneous symmetry breaking in Section IV and the introduction 
of the CTP scheme for classical and quantum systems in Sections V-VI, respectively. The summary of the results is 
presented in Section VII and an Appendix contains the derivation of the CTP Green-function. 

II. SYNOPSIS 

The issue of the time arrow is addressed in this work in a full system of N + 1 degrees of freedom, obeying 
conservative, time reversal invariant dynamics described by the Lagrangian L. We choose a system coordinate, 
denoted by xq, the remaining TV degrees of freedom making up the environment. One encounters two different 
system equations of motion. The Euler-Lagrange equation for xq derived from the Lagrangian L is a non-autonomous 
differential equation. However once the environment coordinates are eliminated by means of their equations of motion, 
this equation becomes an effective, autonomous integro-differential equation. A well known case where it is important 
to distinguish these two kinds of equations is the radiation reaction force in classical electrodynamics [18], sought in 
vain in the elementary equations of motion. 

We follow in this work the change of the elementary equation into the effective one and identify the point where the 
time arrow is generated. For this end it is advantageous to separate two kinds of boundary conditions. Those needed 
to select a unique solution of the effective system equation of motion, called system boundary conditions which remain 
externally imposed parameters of the solution. The remaining, environment boundary conditions are build into the 
effective system equation of motion and are the dynamical source of irreversibility, a dynamically generated system 
time arrow. A necessary condition for irreversibility is to have infinitely many environment initial conditions. 

A simple harmonic model is introduced and its classical dynamics is studied in Section III Classical harmonic 
systems. The model, defined in Section III A System of oscillators consists of linearly coupled harmonic 
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oscillators [19-23] and is describe by the help of the Lagrangian 
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L = 17 X Q ^~ x + J x + 2^ y X n 5— X n - 9n x n x I , (1) 



2 u 2 " w ^ 2 



where a time dependent external source j is coupled to the system coordinate xq. The system trajectory x (t), subject 
to the initial conditions xo(0) = Xo(0) — is constructed in this Section in terms of the external source j(t) and the 
retarded Green-function of the system. 

It is important to distinguish two kinds of time arrow for each degree of freedom. If we impose initial or final 
conditions on a coordinate then the corresponding degree of freedom has a forward or a backward pointing internal 
time arrow. When boundary conditions are used in time, like for variational problems, then no internal time arrow 
exists. Either initial or final conditions are allowed in the model (1) and the retarded Green-function for the system 
coordinate xq {t) contains in its denominator a self energy term, written in terms of the retarded or advanced Green- 
functions of the environment oscillators with forward or backward oriented internal time arrow, respectively. The real, 
dynamical time arrow, called simply time arrow from now on is defined for each degrees of freedom by the direction 
in time in which a linearly coupled source generates a response. In case of independent oscillators the internal and 
dynamical time arrows are identical. But the coupling between the coordinates may change this situation and flip 
or destroy certain time arrows. For instance when an oscillator is coupled to coordinates with different internal time 
arrows its (dynamical) time arrow may simply not exist and the oscillator becomes acausal. Our main interest is 
to understand the rules the system time arrow is defined by the retarded system Green-function, in particular to 
find the conditions for the (dynamical) time arrow be set independently of the internal time arrow, the case called 
irreversibility. Acausality, the nonexistence of dynamical time arrow is a special case of irreversibility. 

We return to a single oscillator in Section III B Time arrow and the null-space where its retarded Green- 
function is written as D r = D n + where D n (—t) = D n (t) and Df(—t) = —D f (t) are called near and far 
Green- functions. It is easy to see that D n , a particular solution of the inhomogeneous equation of motion produces 
the response of the oscillator to an external source in a time reversal invariant manner. Hence the contribution to the 
dynamical time arrow comes from which handles solutions of the homogeneous equation of motion, modes lying 
in the null-space of the equation of motion operator. The lesson is that we have to trace the way the far environment 
Green-functions appear in the solution in our search for the origin of the system time arrow. 

When the internal time arrows are parallel or opposite for two oscillators then the time arrow is well defined or 
nonexistent for the coupled oscillators, respectively according to the results of Section IIIC Coupling of time 
arrows. The state of affairs changes for infinitely many oscillators with continuous spectrum, considered in Section 
III D Continuous spectrum where it is found that even for parallel time arrow the system may not have a well 
defined time arrow. In this case the response of the system for an external source is the sum of retarded and advanced 
terms and the usual causal structure is broken. This happens due to the emergence of poles of the retarded system 
Green-function at complex frequencies on the non-physical half-plane. The imaginary part of the poles yields a finite 
life-time for normal modes, ie. irreversibility. 

A more detailed, dynamical picture, offered in Section IV Dynamical breakdown of reversibility and causal 
structure is based on an analogy with phase transitions. To motivate this view let us start with a simplified version of 
ferromagnetic transition, described by a scalar field (j>(t, x) which is subject to the potential V(<f>) = j<j>+r<f) 2 /2+g(f> 4: /4 
with g > 0. The simple mean-field solution, obtained by minimizing the potential predicts the spontaneous breakdown 
of the discrete symmetry <p — > —(f) for j = as a second order phase transition at r — 0. In the broken symmetry 
phase with r < we find a first order phase transition at j = if the external source j is varied. In fact, the system 
"remembers" the sign of the symmetry breaking term j even after the limit j — > 0. Thought the symmetry which is 
broken spontaneously is discrete this first order phase transition is driven by the soft Goldstone modes arising from 
the breakdown of the continuous symmetry group of the three-space by the presence of domains [24] . 

Phase transitions or spontaneous symmetry breakings take place in infinite systems only and are described in Sec- 
tion IV A Spontaneous symmetry breaking for a realistic, large but finite system as the unexpected emergence 
of a slow collective coordinate. In case of the spontaneous symmetry breaking the space-average of <p(t,x), the order 
parameter displays flip-flops between the minima of the potential in a large but finite system with r < and the 
long time average is the symmetric value (j) = 0. But the typical flip-flop time diverges with the volume and the 
strictly infinite system develops a non- vanishing order parameter. In case of a phase transition without spontaneous 
symmetry breaking the role of order parameter can be played by the proportion of the volume of domains of a given 
phase and the slowing down is experienced without flip-flops which characterize the breakdown of the discrete Zi 
symmetry only. Thus phase transitions are characterized by the non-commutativity of the thermodynamical, N — > oo 
and the long observation time, T — > oo limits. The mathematics is simpler when the averages are calculated as 
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called NT limit, because we deal with a finite system and an unlimited precision of observation. The other way of 
defining averages, 

lim lim N (T), (3) 

T— >oc TV— >oo 

the TN limit, may produce surprising results because certain safely looking arguments of the NT procedure fail due to 
the complications of observing a strictly infinite system within a finite amount of time. Furthermore, in the NT limit 
all information about the system is recorded since there is enough time to resolve the dynamics of a finite system. We 
have limited access to the dynamics in the TN procedure where finite time is allowed to observe an infinite system. 

It is important to realize that phase transition and spontaneous symmetry breaking are the results of such a 
restricted way of observation. The expectation value of the order parameter is vanishing in the ground state or in a 
statistical description of a symmetrical system, a result which can be proven easily in the finite case, N < oo. But 
one should not forget that measurements are always performed in finite time period, T < oo and the large value of TV 
may generate an unexpectedly slow collective mode. For as large N as Avogadro's number some collective coordinate 
may appear to be static for all practical purposes in a Universe of finite age. If this happens with an order parameter 
then the symmetry breaking adiabatic approximation applies with excellent precision. 

The dynamical breakdown of the time reversal invariance is addressed in Section IV B Irreversibility. An 
obvious candidate for the order parameter of irreversibility is the entropy production, S. The definition of the entropy 
is based on limited observability, the separation of the full system into an observed subsystem and its environment 
where only the former can be irreversible. The time scale of slowing down of the order parameter is the time needed 
to restore symmetry by "forgetting" the initial conditions. This takes place when the entropy becomes saturated 
by reaching the heat death, S = and the environment time arrow, set by the intial conditions decouples from 
the system. The formal analogue of the symmetry breaking external source j for irreversible systems is the ie term 
in the retarded Green-functions because the solution of the equation of motion differ by a finite amount when the 
Green-functions with e = + or e = 0~ arc used. But we keep the discussion on a more elementary mechanical level 
and avoid the use of entropy below. 

To exploit the analogy with spontaneous symmetry breaking it is more illuminating to follow the dependence of 
the solution on the order the thermodynamical and the long observation time limits are carried out. Observations of 
finite time can not resolve all spectral levels when the spectrum has a condensation point at the lowest system eigen- 
frequency. As a result, extrapolations based on such limited observations indicate an energy sink in the environment, 
the dynamical basis of dissipative forces and irreversibility. The necessary condition to generate irreversibility in such 
a manner is the presence of a condensation point in the environment energy spectrum at a normal frequency of the 
observed system. One can find similar argument to explain acausality, as well. The time of the action of an external 
perturbation is reconstructed in Section IV C Acausality from observations limited in time and acausality is found 
for sufficiently strong slow soft environment modes. 

Such a heuristic approach to irreversibility and acausality calls for a more systematical description of the effective 
dynamics, constructed for finite time observations. The CTP formalism, used in Sections V-VI is well suited to 
achieve this goal. In the CTP method of classical physics, introduced in Section V Classical closed time path 
formalism the trajectory x(t) for < t < T, subject to the initial conditions x(0) = x(0) — is extended by flipping 
the time arrow at t = T and returning the full system to its initial state for t = 2T, thereby constructing a periodic 
CTP trajectory x CT p(t) = x(t), < t < T and x C Tp{t) = x(T - t) for T < t < 2T. It is easier to handle the CTP 
trajectory by splitting it into two paths, 



X+(t)\ f X C Tp{t) 

x-(t)J \x CT p(2T-t) 



(4) 



The distinguishing feature of the CTP formalism is the independent handling of the trajectories x + (t) and x~(t) that 
might be interpreted as a reduplication of the degrees of freedom. This enlarges the phase space but the equation of 
motion closes this gap by imposing the same dynamics for the CTP doublets. Such a redundancy can be stated in a 
simple manner as a gauge invariance, 

1 + nit) , , . 1 — n(t) 
x(t) = ~^- L x + (t) + -^x- (t) (5) 

reflected by the solution of the equation of motion with an arbitrary function n(t). This scheme admits a variational 
principle for initial condition problems and offers a simple definition of retarded Green-functions even for acausal 
systems. The reduplication of the degrees of freedom makes even possible to find dissipative Euler-Lagrange equation. 
Finally, it separates the oriented and non-oriented interactions among degrees of freedom in time thereby offering a 
clearer insight into the way the environment time arrows influence the system. The action whose variational equation 
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identifies the CTP trajectory is introduced in Section V A Action. The opposite orientation of time for t < T and 
t > T suggests 

Sctp[x + ,x-} = S[x+}-S[x-}, (6) 

where S[x] denotes the action of the full system. The gauge invariance (5) should warn us that this action has a 
large degree of degeneracy. In fact, the action (6) is vanishing for x + (t) = x~(t). The simplest way achieving this, 
together with the regularization of the Green-function with continuous frequency spectrum, is the introduction of an 
imaginary term in the action similar to Feynman's ie prescription. 

The CTP formalism developed in Section V B System of harmonic oscillators solves a riddle about acausality, 
as well. Acausality poses a challenge because the simple integration of the equations of motion in time with an external 
source can not produce an acausal result. But the integration of the equations of motion becomes ill-defined for infinite 
systems due to the possible non-commutativity of the limits At — > and N — > oo. The way out from this ambiguity 
is the use of the Green-function formalism which provides a well defined scheme to handle infinite systems. The 
order of limits (3), to be followed in this case corresponds to periodic trajectories with period length 2T in time. 
Introducing an external source j(t) = jo5(t — t') coupled to the coordinate x generates a causal response in the CTP 
trajectory xcTp(t) for t' < t < 2T — t' . Thus causality appears as a cancellation of the response for < t < t' 
and 2T — t' < t < 2T, a highly non-trivial phenomenon for periodic trajectories. This can be lost in the effective 
dynamics constructed by limited observations where the flipping of the time arrow at t — T can not be performed for 
the unresolved slow environment modes. These continue their motion for t > T as far as the effective system equation 
of motion is concerned. However, the resolved part of the whole system experiences the flipped time arrow for t > T 
and the cancellation needed for causality is destroyed. 

The quantum version of the model (1) is treated within the quantum CTP scheme in Section VI Quantum 
systems, and the same Green-functions are recovered as in the classical case. The formalism is introduced for quantum 
systems in Section VI A CTP formalism and the ways to recognize the classical limit, quantum fluctuations and 
decoherence are pointed out. The issue of the orders of the limits (2), (3) is taken up again in Section VI B Explicit 
breakdown of time reversal invariance because the two orders yield different expectation values even for reversible 
models. It is Fubini's theorem, applied to Feynman graphs which guarantees the same result for both orders (2), (3). 
But a necessary condition for this theorem to hold, the absolute convergence of the loop-integrals is violated by the 
UV divergences of quantum field theories, the natural framework for large quantum systems. One usually follows 
the order (2) but the redefinition of the counterterms, the bare parameters of the theory make possible to implement 
the scheme (3), too. When aiming at detecting a possible irreversibility as a phase transition one should not change 
counterterms at the phase boundary. Rather, one should use a scheme which is applicable to both phases. The 
strategy to save the scheme (2) is based on another way of identifying spontaneously broken symmetries instead of 
following a dynamical procedure, the checking of the slowing down of the order parameter, one carries out a test 
in equilibrium by introducing an explicitly symmetry breaking term in the action and studies the limit when the 
strength of this term approaches zero. A scheme which can be used with the usual procedure (2) for irreversible 
systems consists of performing the limit e — >• in the Green-functions with continuous spectrum, corresponding to 
infinitely long observational time. In fact, a small but finite e generates a finite life-time and represents a dynamical 
realization of the restriction to long time observations mentioned in Section IV B. The small but finite e scheme agrees 
with the usually applied infinitesimal e prescription in the path integral representation. But a non-Hermitean term 
introduced in the Hamiltonian generates non-unitary time evolution and breaks the time translation invariance of 
Green-functions. A simple possibility is proposed which suppresses the non-translation invariant part of the Green- 
functions, the introduction of a degenerate anti-Hermitian operator in the Hamiltonian which generates irreversibility 
by means of friction forces. 

As mentioned before a characteristic feature of the CTP formalism is that it extends the dynamics over a non- 
physical regime by the reduplication of the degrees of freedom. This extension is not cxplorable by the classical 
trajectory due to the gauge invariance (5). Quantum effects, however render them different since x + (t) — x~ (t) 
may be considered within the path integral representation as the source of quantum fluctuations. On the level of 
expectation values, corresponding to the solution of the equation of motion in classical mechanics, one always recovers 
the equivalence of the two doublers, (x + (t)) = (x~ (t)) but one can gain an interesting insight into the effects of 
quantum fluctuations by looking into the interference pattern between the time axes as they build up the observed 
expectation values. Two such effects, related to unitarity of the time evolution and decoherence are discussed in 
Section VI C Unitarity and decoherence. It is shown that the system remains reversible and causal as long as its 
time evolution which may involve memory effects is unitary. Hence the introduction of an infinitesimal imaginary part 
in the Hamiltonian like the one mentioned above is a necessary modification to recover acausality. Another, special 
feature of quantum fluctuations is decoherence, the recovery of classical probabilities. It is shown that irreversibility 
and decoherence are generated by the same dynamical mechanism, both being due to the same spectral strength. 
Their only difference, discovered in case of mixed initial and final conditions for the environment is that the soft 
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quantum fluctuations are weighted with their time arrow for irreversibility and without time arrow for the classical 
limit. 



III. CLASSICAL HARMONIC SYSTEMS 



A simple toy model, a set of coupled harmonic oscillator is introduced in this section with special attention paid to 
the time arrow of each oscillator. The first question is the manner the time arrow is sent through the trajectory of a 
single driven harmonic oscillator from its initial condition. After that the problem of coexistence and competition of 
different time arrows are considered in the simplest case when two oscillators are coupled. Finally the competition of 
time arrows is followed for infinitely many oscillators with continuous spectrum. 



A. System of oscillators 

The Hamiltonian of the model (1), written as 

\ n 71 J n 

shows that this model has a stable dynamics as long as 



(7) 



yJ&-< muj 2 . (8) 

n n 

The solution of the equations of motion 

N 

j = mx a + mulxo + ^ g n x n , 

71=1 

= mx n +mu>lx n + g n x . (9) 

for the time interval < t < T is rendered unique by suitable boundary conditions in time. The degrees of freedom in 
the environment are subject of initial or final conditions, Xj (0) = ij (0) = or Xj (T) = ij (T) = with j — 0, 1, . . . , N, 
respectively. The boundary conditions in time define time arrow for the decoupled system, g n — 0. The coordinate 
Xj with initial or final boundary condition is said to have the time arrow Tj = 1 or tj = — 1, respectively. We use 
initial conditions for the system, tq = 1. 

We seek the system dynamics. For this end the environment coordinates are expressed by means of the second 
equation in (9) as 

SnM = —r, , 9n s 2 aT^oM, (10) 

m[(uj + itT n y - UJ^\ 

where the imaginary term in the denominator with e = + takes care of the orientation of time. This expression is 
inserted into the first equation, resulting the system trajectory 

x(t)= f T dt'D r (t,t')j(t'), (11) 
Jo 

given by the retarded Green-function 

/OO 7 
_ e -Mt-t') D r^ (12) 
-oo 27T 



where 



D r {u) = — fT r-ro 5i (!3) 

v ' m{(uj + ie) 2 - ul\ - £ r (w) v ' 
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with the self-energy 



N 9 



E» = V^ ? : ^. (14) 

n=l 

The roots or the inverse Green-function, Z) r_1 (cj) = belong to the normal mode frequency spectrum. 

When the spectrum of the model N — > oo becomes continuous or at least develops condensation points then it is 
advantageous to use the spectral functions defined for both time arrows separately, 

N 2 

p T (tt) = V 6(tt„)-^K - 0), (15) 

ZTflL)J n 

n— 1 



satisfying the condition (8), 



/ 



fWfl)+^(fl)l<^- (16) 



Models with spectral function p(ft) = O (£l p ) are called Ohmic for p — 1 and dynamics with p < 1 or p > 1 arc 
called sub- or super-Ohmic, respectively [25]. A simple phenomenological form for an Ohmic system with a smooth 
Debey-cutoff is 



pAn) = e(si) _ n _,% ; 02 , , (i7) 



with > and g 2 + g 2 < m 2 ^ 2 ^ 2 -,/^. The corresponding self energy, 

E » = E^/ d0 0f, + 02 (a; + ^2_ 2 5 (18) 



can easily be found, 

TO L—J fl D (u> + ItVLd) ' 

together with the Green-function 



£ r (w) p p • (20) 

m[(w + ze) 2 - w 2 ] + ^(jj^j- - jj^j-) 

It is instructive to compare these results with that corresponding to a simpler spectral function with sharp cutoff 
A [27], 

p > m = ® m J?(lL\\ (21) 



We present the formulae for oscillators with forward internal time arrow only for simplicity, where with 2A p+1 g 2 < 
(p + 1)to 2 cj 2 ^+ 2 which gives for few selected values p = —1,0,1 and 3, 

1 



™f, ,2 , ,2\ , iZ-ng 2 ' 
1 

m(^ 2 -^ 2 ) + ^(2A + z™)' 
1 

m( W 2 - lu 2 ) + _^(f + Ac 2 + i™*) 



D r = \ - 2 • (22) 

TO (w 2 - w 2 )--^ ( ln^ +J .) 
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The equation of motion of models with super-Ohmic spectrum is a differential equation with higher order derivatives. 
The equation of motion of sub-Ohmic systems becomes a more involved intcgro-differential equation with memory 
kernel. The border line between the two families, the Ohmic case, displays the usual friction force and represents 
the only causal and dissipative local equation of motion. Another lesson of the comparison of the self-energy of the 
spectral functions with different cutoffs is that the way the high frequency modes are suppressed may influence the 
low frequency behavior of the Green-function. 

It is sometime advantageous to use the normal modes 

JV 
3=0 

where the mixing coefficients satisfy the sum rule J2 a ^ja = 1 an d yield the retarded Green-function 

A 2 

i[(us + if) 2 — u>, 



° r M = ~u ^k~2 2I ■ ( 24 ) 

v ' m\(i.i A- <iA 2 — r,i 2 l v ' 



The normal mode spectral function 



p n (Q)=sign(n)^SK-|Q|)^|^, (25) 

a a 

defined in terms of the normal mode frequencies oo a allows us to write the retarded Green-function as 

_ f«sa 

y ' J cu + ie-n 

/oo 
dn Pn (n)e- iQt . (26) 
-oo 

Note that whenever there is a condensation point in the spectrum {uj n } then the normal mode spectrum {Sl a } displays 
condensation point, as well. 



B. Time arrow and the null-space 

Let us now take a single oscillator, corresponding to the system coordinate x in our model, decoupled from the 
environment, g n = 0. The simplest way to observe the dynamical time arrow locally, in an arbitrarily small time 
interval around an observation time t' is to find the response to an external perturbation j(t) = jo5(t — t'). The role 
of the boundary conditions in time is to define what is kept fixed when the external source is varied. The time arrow 
point in the direction in which the response to the variation of jo is found. When the response is non- vanishing before 
and after the perturbation then no time arrow exists. The initial or final conditions are usually taken into account by 
the choice of an appropriate solution of the homogeneous equation of motion. Hence the time arrow, defined by such 
test is handled by the null-space of the equation of motion operator D^ 1 = —m(df + Wg), consisting of functions of 
the form Xo(t) — x i cos(uj t) + x 2 sin(uj t), the solution of the source- less oscillator problem. The functions of this 
null-space are free fields in classical electrodynamics and mass-shell modes in quantum field theories. We shall call 
modes inside and outside of the null space free and driven modes, respectively. The variational equation of motion 
determines the driven modes only because the free modes drop out from the Lagrangian [26] and are fixed by the 
boundary conditions rather than dynamical equations. In fact, the response of the oscillator to an external source is 
ill defined, it diverges at resonance when the frequency of the source approaches luq. 

The usual way of separating the null-space and the rest is the introduction of far and near Green-functions by 
analogy with electrodynamics, 



D f {uj) = -^sign(cj),5(a; 2 - u 2 ), 



D%{lo) = P I (27) 



where P denotes the principal part. The time-dependent forms are 

2mw 
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The far or near Green-function handle the free and the driven modes, respectively and the retarded and advanced 

r r 

Green-functions are Dq — Dq ± D J . The driven modes are generated by the given external source but have no 
time arrow. In fact, an external source j(t) = jod(t — t') generates a response both for t < t' and t > t' because 
the near Green-function is symmetric, D^(t) = Dq{— t). The time arrow is introduced by the antisymmetric far 
Green-function, D^—t) = —D^t). 

The lesson of these trivial remarks is that we have to separate the components of the null-space from the rest in 
discussing the issue of the time arrow. A time reversal invariant linear equation of motion generates a symmetric 
response for the external source but leaves the issue of the null-space open. The mathematical source of an eventual 
anomaly of the time arrow is a strong singularity of the near Green-function at the null-space which manages to mix 
these two components. Such a mixing takes place when oscillators with different null-spaces are coupled because the 
non-dynamical free modes of one oscillator are coupled to the dynamical, driven modes of the other. This is the 
mechanism by which the environment boundary conditions influence the system time arrow. 



C. Coupling of time arrows 

What happens with the time arrows in our model, defined for g n = when the oscillators become coupled? The 
answer is presented for the simplest case of N = 1, when the retarded system Green- function 

D r (oj) = 5 (29) 

m[(w + it) 2 - c 2 ] - % (u+u * )a _^ 

has four poles, u a a > = auj a i + ir/ a > with a, a' = ±1. Due to the stability condition (8) av and r) a ' are real and 




1 T. — 1 

(30) 

with d = \J (wq — ^i) 2 + 4g 2 /m 2 . In case of parallel decoupled time arrows the interaction does not change the 
time arrow as expected. But the environment destroys the system time arrow when the decoupled time arrows are 
conflicting. When the two null-spaces agree, ojq = w\ and the time arrows are conflicting then the poles remain on 
the real axis and the conventional Green-function method fails. 

The generalization of these results for larger, finite N can easily be imagined though the details are rather involved. 
The null-space of the normal modes is 2N dimensional and these free modes are fixed by the boundary conditions. The 
contribution of the normal modes can easily be seen by using the partial fraction decomposition for the Green-function 
(13) where each partial fraction corresponds to a normal mode. If all degrees of freedom have forward decoupled time 
arrows then the stability condition (8) makes all poles of Dq(uj — ie) real and there is a global time arrow. In case of 
conflicting time arrows in the decoupled case the normal modes possess well defined orientation in time in the absence 
of real poles. But some accidental degeneracy of the decoupled null-spaces with opposite time arrow may create ill 
posed problem. 

What happens in the limit N — > oo? There should be no qualitative change as long as the spectrum remains discrete, 
without condensation point. The distinguishing feature of the discrete spectrum is that the analytical continuation 
of Green-functions over the complex frequency space is unique. This property, together with the stability condition 
(8) keep the normal frequency spectrum real and the observations, made for finite systems apply. But a condensation 
point or a continuous part in the spectrum may lead two important changes. When environmental free modes with 
opposite time arrow converge in frequency towards a system free mode then for strong enough environmental spectral 
strength the minority system time arrow may simply be erased. This phenomenon, irreversibility, is the subject of 
the rest of this work. 
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D. Continuous spectrum 

When the environment spectrum has a continuous spectrum then the spectral functions (15) yield the retarded 
system Green-function 

^ U i • \9 91 9+^D- (Or)_+iw)+g^O D+ (n D+ -iw) i^) 

m\(LO + ie) z - ufA + tt + n — „ ? — — n n — 

/ UJ mUD+SiD- (w+ii2r>+)(w— i!2d- ) 

The poles of the self energy, given by Eq. (14) merge and accumulate their imaginary parts in Eq. (19). This self 
energy generates poles in Eq. (31) with finite imaginary part when inserted in the Green- function (13), indicating 
that normal modes are damped and have finite life-time, the microscopic manifestation of irreversibility. Note that 
the Green-function is in general acausal, there are poles on the unphysical sheet. Poles with finite imaginary part, 
leading to reversibility are present for spectral weights with different suppression for large frequencies suggesting that 
the breakdown of the time reversal invariance arises from the behavior of the spectral strength at vanishing frequency 
and reversibility may prevail for spectra with a single condensation point, too. 

In the case of a pure initial condition problem, = the local time arrows, defined for g n — agree. Nevertheless 
the Green-function (31) may be acausal, the condensation point in the spectrum of forward moving modes can generate 
a backward propagating signal. We shall return to the origin of this rather surprising effect later. 



IV. DYNAMICAL BREAKDOWN OF REVERSIBILITY AND CAUSAL STRUCTURE 

We set identical, forward time arrow for each oscillator from now on and turn to the question of the breakdown of 
the time reversal invariance of the effective system dynamics in the thermodynamical limit N — > oo. For finite N the 
effective system equation of motion, obtained by eliminating the environment coordinates is a differential equation of 
order 2N, needs 2N boundary conditions and remains reversible. The effective equation of motion continues to be 
reversible in the limit N — > oo for discrete spectrum due to the availability of the normal mode decomposition and the 
Green- function (26). But this state of affairs may change when the spectrum has a condensation point and infinitely 
many poles coalesce as in the Green- function (31) with g_ = 0. The effective equation of motion needs less boundary 
conditions then expected and the remaining "invisible" boundary conditions arc the source of irreversibility. It is 
argued below that this mechanism is similar to spontaneous symmetry breaking. 



A. Spontaneous symmetry breaking 

The Green-function (13)-(14) is unique with well defined analytical continuation over the complex frequency space 
for discrete spectrum and the order of the summation in the self energy (14) and the integration over the frequency in 
Fourier integrals to find the time dependence of the Green-function is arbitrary. But if there is a condensation point 
in the spectrum then these operations may not commute anymore. The relevance of the order of the limits where 
the system size and observation time is sent to infinite is a hallmark of spontaneous symmetry breaking. Note that 
the normal mode Green-function (26) remains always causal because the identification of the normal modes always 
precedes the frequency integration. 

What is the correct order of limits? Mathematical ambiguities always point to a choice in the preparation or in 
the observation of the system and the choice between the NT and the TN limits, given by Eqs. (2)-(3), respectively 
depends on the physical circumstances. It is instructive to consider the status of rotational symmetry in a macroscopic 
body of characteristic length scale I = 1cm, being composed of N ~ 10 24 atoms of mass m whose elementary 
interactions are assumed to be invariant under translations and rotations. We separate the translational and rotational 
motion by distinguishing laboratory and body-fixed, co-moving coordinate systems, the latter is defined by having 
vanishing center of mass velocity and diagonal tensor of inertia 

»> L ( 32 ) 

n 

Denote the translation and rotation which bring the laboratory frame into the body-fixed coordinate system by T(X) 
and R{6, <p, a), respectively and use the coordinates X, the Euler angles 0, <p and a and 3iV — 6 relative coordinates 
expressed in the body-fixed frame to describe the positions of the particles of the body. The effective dynamics of the 
three Euler angles, considered as order parameters for rotation is driven by the Hamiltonian 



H rot = ±Li(Q-y k L k + 0(L'), 



(33) 
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written in terms of the total angular momentum L. 

One can show that the ground state of a finite system described by a single component wave function and obeying 
the Schrodinger equation with regular potential is non-degenerate [28] . This result makes the angular momentum in 
the ground state of arbitrary large but finite bodies, held together by rotational invariant forces vanishing and excludes 
spontaneous symmetry breaking. The resolution of this apparent contradiction is provided by the adiabatic, classical 
approximation which breaks rotational symmetry for N — > oo and is an excellent approximation for as large N as 
Avogadro's number. The Euler angles of a finite body follow classical dynamics in a good approximation because the 
large inertia tensor makes the rotational excitation spectrum dense. The low lying rotational excitations decouple from 
the remaining, internal degrees of freedom and the ground state can be constructed by means of angular momentum 
eigenstates whose number is approximately N- independent giving (Lj) = hO (AT ). The typical angular velocity is 
therefore of the order of magnitude 

where the mass of the proton is used for m. The exceedingly large time for a complete turnaround makes the rotation 
unobservable and the macroscopic body appears to break rotational symmetry in its ground state. 

Note the important role intial conditions play in this mechanism. Spontaneous symmetry breaking takes place when 
the ground state or the ensemble average display less symmetry than the underlying equations of motion, when the 
initial conditions prevent the system to explore the phase space prescribed by symmetry. If more energy is available 
in the initial state then the order parameter may move faster and the symmetry may appear to be restored. In other 
words, critical slowing down prevents the system to "forget" its initial conditions which were imposed at as low energy 
as possible. 



B. Irreversibility 

Our goal is now to understand irreversibility as a natural result of the widening of the discrete spectrum lines due 
to insufficient observation time [27]. Let us return to the question of the issue of the limits NT and TN. Observations 
made according to the NT procedure can resolve every excitation frequency, can trace the energy exchange among 
coordinates and can recognize the detailed, microscopically reversible dynamics. This is not the case for the TN 
procedure where no finite observation time can resolve all frequencies around a condensation point. To see this point 
clearer let us follow the system for a finite amount of time and look into its response to the source j(t) = joS(t — t'). 
The finiteness of the observation time is taken into account by considering the product 

x° bs (t) =c(t)x (t) (35) 

as the measured result where c(f + t) = c(t' - 1) is an IR cutoff, c(t' + 1) = 1 for \t\ < T and c(t' + 1) = for \tT > T. 
The observed frequency spectrum 

/OO 7 ' 
— c(w - lo'W\J) 3 {J) (36) 

can be written as 

P ° bs (n) = V-^[ c (n-u a )-c( Wo + o)] (38) 

* — ' 4irmu! a 

a 

denotes the spectral density extracted from the limited observations. If there is time to resolve all spectral levels, 
c(uj a — cjp) « a ^ (3 then we find peaks, 

p"bs {n) „ signin) J2^-c(n-u; a ) (39) 



where 



whose separation is more obvious by making observations in longer time. But there is always an unresolved part of 
the spectrum with condensation point and the corresponding infinitely many normal modes provides the energy sink 
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needed for dissipation. An IR cutoff on observations introduces coarse graining and that generates irreversibility [29] . 
The spectral lines spread even if the source j(t) = jo8{t — t") with t" =/= t' is located in an asymmetric manner in time 
with respect to the IR cutoff except that the speed of spreading depends on t" — t'. 

The characteristic feature of irreversibility, generated by the spontaneous breakdown of the time reversal invariance 
is that the action is real and poles of the Green-functions with non-vanishing real part appear in pairs, w± = ±cji+zcj 2 - 
The real part, ±cji is the remnant of the formal time reversal invariance of the original dynamics and the imaginary 
part, iuj 2 reflects irreversibility. 



C. Acausality 



It was pointed out above that the Green-function (31) is acausal even if g- = when initial conditions are imposed 
for all degrees of freedom. The origin of this surprising effect can also be understood by taking into account the 
finiteness of observation time. We use now the observed trajectory (35) to reconstruct the source 



j° bs (u) = D r - 1 (cj)x° bs (Lo) 



(40) 



written as 



b f dio , D r (uj ) , 

j obs (oj) = / c(w-u ) -7-4- 3 w). 

J v ; J 2tt v ' D r {uj) Jy ' 



(41) 



The partial fraction decomposition can be used to obtain a Green-function (13)-(14) as a sum of simple pole terms, 



D r {uj) = y __£p_. 



(42) 



The contribution of the term p dominates the integral in Eq. (41) for lj' ~ co p , 



dw' 



IV' — bJ. D 



(43) 



We choose the IR cutoff c{u) = 2j]/(uj 2 + rj 2 ) with rj = 2n/T and find for the source j(t) = jo5(t) 



jT(t)=h 



S(t) - ^ e ~i^-^ 



(44) 



Apparent acausality results from the difficulties in reconstructing a sharply localized source in time due to the abundant 



soft modes with w p w 0. 



V. CLASSICAL CLOSED TIME PATH FORMALISM 



It was pointed out in Section III B that the driven modes build up an unoriented motion in time and the time arrow 
is set by the free modes. It is easy to recast the equation of motion in a slightly modified form where these two classes 
of modes are clearly separated. Apart of making this point explicitly, this formalism is needed, as well, to tackle 
difficulties of infinite systems. The point is that the solution of the equations of motion of a finite system can easily 
be obtained by direct integration in time for N < oo but complications arise when spontaneous symmetry breaking 
is suspected for N — > oo because both exact solutions and the numerical quadratures belong to the NT scheme. A 
simple formalism presented in this Section allows us to treat infinite systems by means of powerful functional methods 
(Green-functions) in either the NT or in the TN scheme. 

One needs the inverse of differential operators when a degree of freedom is eliminated by its equation of motion. It 
is a well known problem that the self-adjoint extension of the derivative operator id/dt requires well defined boundary 
conditions, not available in an initial condition problem. How to save the functional method in this case? The existence 
of Green-functions follows from the functional setting of the action principle. In fact, let us take a dynamical system 
described by the coordinate x and the action, 



S[x] = ^xD^x + O (x 3 ) , 



(45) 
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written in condensed notation, introduce an infinitesimal external source j(t) for diagnostic goal and define the 
functional 

W\j]=S[x]+jx (46) 

in which the trajectory x is chosen to satisfy the variational equation of motion, 

6S[x ] 
5x 

The form 



+ j = 0. (47) 



W\j] = ~jDj + 0(f), (48) 

identifies the Green-function as a second functional derivative, 

S 2 W\i] 

D = ?-£M . (49) 

This procedure can not handle initial conditions in time because the quadratic form D _1 of the action must be 
symmetric, D — D n , excluding the null-space modes in the solution of the equation of motion. Another face of this 
problem is that the variational equation must be imposed at the end point t = T if an initial condition problem is 
considered and this condition cancels the final momentum. To avoid this restriction we flip the time arrow at the end 
of the motion and return the system backward in time to its initial state. It turns out that such an extension of the 
motion leads to non-symmetrical Green- functions, such as D r in the variational equations. The trajectory xcTp(t), 
< t < IT obtained in this manner can be broken into a pair of trajectories with different time orientations [26] as 
in Eq. (4) and this formalism goes over to the CTP method [6] when considered in quantum mechanics. 



A. Action 



There are two different realizations of the CTP action principle, one for T < oo and another for T = oo. The limit 
of T — > oo of the first is not continuous, it does not produce the scheme T = oo because the former has no null-space 
in the quadratic form of the linear equation of motion. Let us start with the scheme T < oo by introducing two 
trajectories for each degree of freedom, x(t) — > (x + (t), x~ (t)) = x(t), satisfying the boundary conditions x ± (0) = 0, 
i ± (0) = Vi, x + (T) = x~{T) and the classical CTP action is defined as 

Sctp[x]= ( dt[LcTp(x + (t),x+(t))-L* CTP (x-(t),x-(t))}, (50) 
Jo 

#tLcTP = L being the original Lagrangian of the model. The minus sign in the action reflects the opposite orientation 
of time for the two CTP trajectories. When the original, real Lagrangian is used then the CTP action is vanishing 
whenever x + (t) = x~(t). This degeneracy is due to the lack of initial condition for the velocity. Such an initial 
condition could in principle be imposed by using the extended action S[x] — > S[x] + Sb c [x] with 

\ ox |t=o ox 

where £ is a Lagrange multiplier and pi denotes the initial generalized momentum. A simplified scheme is used below 
with Vi = corresponding to I = for theories with space inversion invariance. The degeneracy of the CTP action for 
real Lagrangian and I = for x + (t) = xr (t) can easily be lifted by making a small imaginary shift of the harmonic 
potential, 

L C tp(x,x) = L(x,x) + i^x 2 (52) 

for a particle of mass m and the limit e — > is performed after deriving and solving the variational equations of 
motion. One could have used real shift of the potential but the the imaginary shift is closer to the way the degeneracy 
is lifted for T = oo. 
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Let us consider the simplest example, a harmonic oscillator in the presence of a driving source j when the CTP 
action can be written as 



(53) 



Sctp[x] = J dt 

with j = (j, —j). The equation of motion x = —Dj for the physical trajectory is 

• T W = - E C dt ' Daa ' (*' O^O. (54) 
a , Jo 

for < t < T where D aa = (D) aa denotes the matrix elements of the CTP Green- function with a, a' = ±1 and a is 
arbitrary. The independence of the trajectory (54) of the choice of a yields the CTP relation 

D++ + D—=D~ + +D+-, (55) 

assuring that the Green-function can be parameterized by three functions, 

h (D + D n D-Df\ . . 

D ={D + Df D-D-)- (56) 

The inverse Green-function D^ 1 is a symmetric operator, D~ ltr = D^ 1 hence so is D and D(t,t') = D(t',t), 
D n {t,t') = D n (t',t) and D f (t,t') = (*',*). According to Eq. (54) the retarded Green-function is D r = G n + G f 
and the advanced Green-function will be defined by D a = D n — D* . The function D(t,t') can not be identified by 
the help of Eq. (54). The explicit calculation of D n and is given in Appendix A3. 

The limit T — > oo of the CTP Green-function is discontinuous because the coupling of the two time axis, introduced 
at t = T, disappears if T = oo and the Green-function recovers translation invariance in time. How to maintain 
the coupling between the two time axis, the distinguished feature of the CTP scheme? This coupling is supposed to 
handle the null-space modes of the linearized equation of motion hence one expects the form 



/oo />oo 
dtL(x + ,x+(t)) - / dtL(x-,±-(t)) + S B c[x} 
-OO J —OO 



(57) 



where a translation invariant Sbc[%] = O (x 2 ) singles out the null- modes which couple the two time axis. The 
determination of the Green-function starts with the observation that the Fourier transform of D n and D^ defined by 
Eq. (56) are real and imaginary, respectively. This property, together with the identity 

' P- - inS{x) (58) 



x + le x 

requires the usual form, Eqs. (27), for D n (uj) and D^(uj). A simple form of D which reproduces the desired D n and 
D* is 

% / e -iu>o\t— 1'\ e -iuo(t-t')\ 

= Ve <Wo(t "° e^l t+t 'l ) ' (59) 

and a formal prescription to derive it as T — > oo is presented in Appendix A 4. We are interested in a trajectory 
for t ~ T/2, safely away from the endpoints t = and t — T which should decouple from the dynamics to recover 
translation invariance in time as T — >• oo. In order this to happen we would need short enough life-time, wo/e C T, 
generated by the imaginary term in the action for the pole contributions. However such a removal of the dependence 
on the end points leads to decoupling of the two time axis taking place at t = T and thus the CTP scheme is lost. 
Therefore we have to rely on the cancellation due to the fast rotating phase of the pole contributions as T — > oo to 
recover translation invariance in time. The usual way to realize this scheme is to treat e as a formal, infinitesimal 
quantity. 

The discontinuity of the Green- functions in reaching the continuous spectrum inherent in setting T = oo, is 
reminiscent of phase transitions and in this respect the continuous spectrum and T correspond to the adiabatic 
approximation and the system size, respectively. The realistic observations belong to the scheme of T < oo, containing 
the scale parameter T. If T is well beyond the observation scales then it is more natural to work out formal rules for 
the case of continuous spectrum where this scale parameter is absent, and to use them to describe observations made 
at large enough T. 
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The inverse of the Green-function, calculated in Appendix A 5, 



D 1 (uj) = m 



(- 2 --o 2 ) 



it 



1 

-29(w) 



-20(-u>) 



(60) 



is the quadratic part of the translation invariant CTP action for T = oo. It shows that the coupling between the time 
axis interchanges modes with given time orientation and the time arrow of the solution is set by the "interference" 
between the time axes. The inverse Fourier transform of this expression yields 



Sbc[x] = -P f 

J-< 



dtdt 



,x + (t)x-(t') 

t-t' 



f 

J — C 



dt[x+{t) -x-{t)] 2 



(61) 



where the recognition of the time orientation of the modes requires a non-local term. The coupling of finite strength 
between the time axis at t — T of the scheme T < oo is traded into an infinitesimally strong coupling in the latter, 
smeared out over an infinitely long time evolution. 

We have two different ways of introducing the intial conditions. One possibility is to set up initial conditions at 
some initial time ti ~ T/2 by means of Lagrange multipliers. For such a trajectory the integral on the right hand side 
of Eq. (54) extends between ti and infinity and D n + D' with continuous frequency spectrum is just the retarded 
Green-function as expected. Another possibility is to start the motion with Xi — Xi — at t = and let the source j 
drive the trajectory to the desired initial conditions used before at ti ~ T/2. As the source j varies for < t < ti x(t) 
and x(t) sweep through the set of accessible initial conditions. For such initial conditions the two realizations of the 
motion are equivalent and we can use the simpler scheme, starting the motion at U and using translation invariant 
Green-functions. By enlarging the class of dynamical quantities which enter in the Lagrangian with an external source 
we can always describe the actual experimental preparation of the system and render our intial conditions accessible. 
This is the way the restriction Vi = on the initial condition, mentioned after Eq. (51) can be avoided in the CTP 
formalism. 



B. System of harmonic oscillators 

The steps presented in Section (III A) to find the system trajectory can be repeated within the CTP formalism 
without difficulty. Eqs. (9) then refer to CTP pair of trajectories, the inverse system Green-function is of the form 

D- 1 = D 1 - ata (62) 

where the Dq 1 is given by Eq. (60), the matrix 



1 

-1 



(63) 



represents the different signs in front of the two actions in Eq. (50) in the coupling between the oscillators. The CTP 
matrix of the self-energy, 



SM^^S^H, (64) 



71=1 



contains D n (uj), the Green- function of the n-th environment oscillator which is calculated in Eq. (A28) for n = 0. 
The spectral representation for the self-energy is 



poo 

±(ui) = V / dn2Slp T (Sl)D (TW,Sl), 

r=±l J 



(65) 



D(uj, f2) being defined in Eq. (A28) for T — > oo. The off-diagonal elements in (A28) can be used to find the useful 
relations 

eh = -ttMI^d + MM)], 

E'(w) - -iimga{u)\p + (\u\) - p-(\u\)]. (66) 
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It is an important property of the CTP Green-functions that D, uDq a and £ all have the block structure shown 
in Eq. (56), allowing to define the functions D n , D~ ln , E™, etc. One carries out the spectral integral first in the self 
energy within the scheme (3) and the limit T — > oo follows after that. The spectral weight (17) results 



S "( W ) = --Y 
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2 



£/(„) = -i™ y 



^ - E u^sry 1671 

The inversion in Eq. (62) can easily be carried out [16] and one finds 

D n = D r D- ln D a 
D f = -D r D- lf D a 

D l = -D r D- u D a (68) 

in particular 

D a = — ± -, (69) 

D^ la - n a 

in agreement with Eq. (13) and the results of Section III A are reproduced within the CTP formalism. Despite 
this equivalence the CTP formalism offers advantages when effective equations of motion are dissipative. In fact, 
it is known that variation method can not lead to non-conservative dynamics. But the reduplication of degrees of 
freedom offers a possibility to derive from the action principle, cf. the odd powers of oj in the denominator of the 
Green-function (31). 

Another advantage is that this formalism offers a mechanism explaining the loss of causality for dissipative systems, 
noted in Section IV. The price of transforming an initial condition into a variational problem is to make both the 
trajectories and external perturbations periodic in time. The response on the trajectory xcrp(t), used in Eq. (4) to 
a perturbation made at time t' is restricted to the time interval t' < t < 2T — t' . Thus causality becomes the result 
of a rather fragile cancellation within the time interval 2T — t' < t < IT when the equation of motion is integrated 
from t = to t = 2T and inaccuracies like the one mentioned in Section IV C may lead to acausal dynamics. 

One may get an insight into such a breakdown of causality by looking into the "polarization cloud" the system 
generates in the environment, the correlation between the system and the environment motion. This correlation 
represents the "dressing" of the system particle by the environment, the emergence of a system quasi-particle. When 
the time arrow of the effective system dynamics is flipped at time T in the framework of the CTP formalism then 
the motion of the environment which is related to the system is reversed, as well. But irreversibility, the spontaneous 
breakdown of time reversal invariance implies a limited resolution of the soft slow environmental modes by the system. 
Such a limited resolution leads to an incomplete reversal of the environment motion at t — T within the effective CTP 
dynamics. In other words, that part of the "polarization cloud" which belongs to the very soft modes and remain 
unresolved within the limited observation time continues its motion and what is turned back at t = T is actually 
only the observed fraction of the dressed quasi-particle. It is natural that this deformed quasi-particle follows a path 
backwards in time which is different from the forward motion and the cancellation for 2T — t' < t < 2T is incomplete. 



VI. QUANTUM SYSTEMS 

The system of quantum harmonic oscillators, defined by the Hamiltonian (7) can be solved in a manner analogous 
to the discussion of Section III A. The equations of motion (9) are imposed on the coordinate operators taken in 
the Heisenberg representation and the solutions can be written by means of the same retarded Green-function as 
in the classical case leaving the result of Sections III-IV valid for quantum oscillators. Nevertheless one gains more 
insight and efficiency for interactive, anharmonic models when the CTP formalism is followed in the path integral 
representation. 
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A. CTP formalism 



The quantum mechanical expectation value 

im\A\m) = y>(o)\u*(t,o)Au(t,o)Mo)), (70) 

where \ip(0)) denotes the initial state and U(t, 0) stands for the time evolution operator. It is crucial to realize that 
one has to treat the operators U\t, 0) and U (t, 0) separately in the expectation value to build up eventual interactions 
in perturbation expansion and to generate the observable A by functional derivatives. Such an independent treatment 
is possible if the source j is allowed to be different in the two operators. Such a non-physical, formal extension of the 
formulae of the expectation value not only allows us to derive Green-functions for such an initial condition problem 
but gives access to the density matrix, as well. The basic idea is therefore to deal with two time axes in an explicit 
manner, one for the operator U and another for W. Such a point of view proposed in the original CTP formalism 
[6], based on the generator functional 

e iwB] = Tr[C/(j+;T,0)p(0)[/t(-j-;T,0)]. (71) 

where U(j; T, 0) is the time evolution operator in presence of the external source j coupled linearly to the coordinate 
and p(0) denotes the initial density matrix. The physical system with unitary time evolution corresponds to the choice 
j + = It is not difficult to find the path integral representation for this functional, 

e iw\j] = f D[^ e is TPlx]+ifdt3{t)x(t) } (72) 

where the action is defined by Eqs. (50) and (52) and the integration is over trajectories with fixed initial points, 
x ± (0) — Xi and arbitrary common final point, Xf = x + (T) = x~{T). 

The expectation value of the coordinate is given by the functional derivative, 

The equation of motion for the expectation values belonging to the quantum mechanical intial condition problem, 
can be derived from an action principle, just as in the classical case. One performs a functional Legendre transform 
by defining the effective action, 

T[x] = W[j] - xj (74) 

in condensed notation where 



Sj 

The inverse Legendre transformation is given by (74) and the change of variable 

ST[x] 
3 Sx ' 



(75) 



(76) 



whose defining equation is interpreted as an equation of motion. Note that Eqs. (45)- (47) define similar Legendre 
transformation in classical mechanics, leading back to the classical action as effective action. Naturally the effective 
action is different from the original one when some environmental degrees of freedom are eliminated, and effective 
dynamics can be considered both in classical and quantum cases. 

It is instructive to extend the CTP formalism to an Open Time Path (OTP) scheme, based on the generator 
functional 

e *m?;x±] = {x +\ U{j +.T,0)p(0)UH-j-;T,0)\xJ) (77) 

to be interpreted as the matrix element of the density operator (x~f\p(T; j)\xj), subject to a generalized time evolution 
which has external source j + and j~ for U and U' , respectively. The corresponding path integral expression, (72) 
contains open trajectories, with common initial point, 2^(0) = Xi and independent final points, x ± (T) = x^ . Note 
that the integrand of the path integral is the product of factors, one depending on x + and the other one on x~ if 
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the initial state is pure, the initial density matrix is factorisable and all degrees of freedom arc present. The two 
trajectories of the pair x are uncorrelated in the path integral and the OTP formalism produces in this case simply 
the product of a transition amplitude and its complex conjugate in one expression. 

The formalism with double time axes becomes advantageous when mixed initial state is considered or some degree of 
freedom is eliminated because the intuitive scheme of Feynman diagrams can be maintained to represent expectation 
values in terms of elementary processes. Though this issue is not important when harmonic systems are considered 
there is another advantage which is relevant here. The integrand of the path integral is not factorisable when mixed 
states contribute to the expectation value and this makes the trajectories of the pair correlated within the path 
integral. We can follow the quantum-classical transition by inspecting such correlations generated by the CTP action 

Sctp[x] = S[x+] - S*[x-\ + Si[x], (78) 

which is not additive due to the influence functional [13] Si[x}. The integrand of the usual path integral with a single 
time axis is the contribution of a given trajectory to a transition amplitude [30] . In a similar manner that integrand of 
the OTP path integral represents the contribution of the pair of trajectories, x to the density matrix. Its suppression 
for x^ ^ xj signals decoherence in the coordinate basis [31, 32]. When all degrees of freedom arc present then the 
integrand is a pure phase and its suppression may come from fast rotating phase as the trajectories are varied. When 
degrees of freedom have already been eliminated then SsSi [x] > represents environment induced decoherence. 

The two trajectories of the pair are always identical, x + (t) = x~(t) for the solution of the equation of motions in 
the classical CTP formalism. The solution of the classical equation of motion corresponds to expectation value in 
the quantum case and the time axes give identical time evolution for the expectation value of the coordinate. But as 
soon as we look into the interference pattern between the time axes quantum fluctuations appear in x + (t) — x~ (t) [33] 
witnessing incomplete decoherence and indicating that the possibility of separating the two trajectories of the CTP 
pair is an O (h) genuine quantum feature. 



B. Explicit breakdown of time reversal invariance 

We make a little digression in this Section and consider a technical problem in the procedure of detecting irre- 
versibility, generated by the well known divergences of quantum field theory. The question of the order of the limits 
NT-TN is specially acute for the perturbation series of interactive quantum systems where it corresponds to the 
choice of the order of summation over the energies and the momenta of intermediate states. In fact, Fubini's the- 
orem guarantees the independence from the order of energy and momentum integration for finite Feynman graphs 
only. Divergent graphs change their values when the order of integration is changed. One usually follows the NT 
scheme and integrates over energies first because this scheme allows the proper treatment of infrared divergences 
when necessary. One may use the scheme (3) as soon as the thermodynamical limit is proven to be safe, the price 
being a modification of the counterterms controlling the ultraviolet dynamics. But rather than making the change 
NT -o- TN at a phase transition it is more preferable to rely on the same NT procedure in exploring the whole 
phase diagram. Therefore we can not avoid the problem of establishing spontaneously broken symmetries within this 
scheme. The usual solution, mentioned in Section II is based on equilibrium observables, one breaks the symmetry 
by an external source j, coupled to the order parameter. Since the symmetry is broken explicitly for a fixed j ^ the 
order parameter assumes a constant, nontrivial value and the scheme NT is obviously applicable. The symmetry is 
recovered or remains broken spontaneously after the limit j — > has been performed in the symmetric and the broken 
symmetry phase, respectively. 

The simplest way to implement this procedure for irreversible systems is to keep the ie term in the denominators of 
the Green-functions small but finite rather than infinitesimal. The imaginary part of the poles of the Green-functions 
spreads the discrete spectrum lines and the resulting finite life-time of excitations represents a dynamical realization 
of the IR cutoff for observations, introduced by hand in Section IV B. Now we have three limits to perform, N — > oo, 
T — > oo and e — >• 0. Since the ie term regulates the pole contributions to the Green- functions the limit e — > must be 
carried out after the other two limits. One can define in this manner the limits eNT, 



and eTN, 



lim lim lim O eN (T), (79) 

e-s-0 JV-s-oo T-s-oo 



lim lim lim O eN (T), (80) 



and either of them can be used for irreversible systems. 



19 



It should be noted that neither of the limits eNT or eTN is equivalent with the usual implementation of Feynman's 
ie prescription because e is not infinitesimal. To see the difference between small but finite and infinitesimal e let us 
briefly outline a slight generalization of the procedures eNT or eTN where irreversibility is achieved by introducing an 
O (e) non-Hcrmitean term in the Hamiltonian H = H{t) — iK with H^(t) = H(t), = K where H{t) might be time 
dependent due to some external sources. It is the non-unitarity of the time evolution which introduces complications 
in setting up the Heisenberg representation where the equations of motion are easier to access. Let us start with 
the Schrodinger representation where operators As are time independent and states are subject of the Schrodinger 
equation ihd t \tp(t))s = (H(t) — iK)\ip{t))s- One can define a diffusive Heisenberg representation by performing the 
basis transformation 

A d (t) = U\t,U)A s U(t,U) (81) 

with 

U(t,t') = T[e~^ f n dt ' H{t,) } (82) 
which places the unitary time dependence in the operators and generates the equations of motion 

imm)) d = -iK d {t)\ip(t)) d 

ihd t A d {t) = [A d {t),H d {t)\. (83) 
The expectation values satisfy the equation of motion, 

d t m)\Mt)\m) = ^(m\[Mt),H d (t)]-i{A d (t),K d ( t )}\m) m 

showing an unexpected effect of non-unitary time evolution, the loss of invariance under translation in time for 
Green- functions. This happens because the factors of the time evolution operators corresponding to non-Hermitean 
Hamiltonian in the usual Heisenberg representation do not simplify in the Green-functions according to the usual 
rule U(t n +i,ti)U'(t n ,ti) — U(t n +i,t n ). A more physical way to see this is to note that the introduction of a non- 
Hcrmitcan part in the Hamiltonian is a dynamical way to implement an IR-cutoff like the one used in in Section IV B 
because 1/e is proportional to the life-time of the states. Since all state decay the initial time when the initial state 
with a given norm is set remains always explicitly present in the expectation values. Irreversibility is generated even 
in the scheme eNT as long as e ^ because the integration contour of the Green-functions runs at finite distance 
from the poles what can be interpreted as a "natural", finite width of the spectral lines. 

It is instructive to find the analogue of Ehrenfest's theorem in this representation. For this end we choose the 
Hamiltonian 

H(t) = £- + V(x)-j(t)x, (85) 



2m 



and use K = \he. The equations of motion 



dt 

can be written as a Newton equation with a friction term, 



d i \ P / \ 
— (x) = e(x), 

lip) = (j-V'(x))-e(p), (86) 



dt 2 ' ' IJ x el}y " dt 

with an effective potential 

Veff(x) =V(x) + ^me 2 x\ 



m i( x ) = <K//(z)> - 2em^-(x) (87) 



reflecting a suppression of deviation from the initial position by the decay of the state norm. 
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C. Unitarity and decoherence 

Let us return to our simple harmonic model denned by the Lagrangian (1) where the effective system dynamics is 
described by the the generator functional 

= e-M®, (89) 
of the form identical to the classical counterpart shown in Eq. (48). It gives rise to the Green- function (49), 



^i-(*)^i + (* , )| 3=0 



(0\T[x(t)x(t')}\0) = ihD ++ (t,t'), 
(0\T[x(t')x(t)\0)* = ihD—(t,t'), 
(0\x(t')x(t)0) = ihD + -{t,t'), 

(0\x(t)x(t')0) = ihD~ + (t, t'), (90) 



defined by following the convention of including some ih factors. Due to the identity of the functional W[j] for 
classical and quantum harmonic systems, Eqs. (48) and (89), respectively, the classical and quantum CTP Green- 
functions agree. The Legendre transformation (74)-(75) leads back the the classical action for harmonic system and 
the argument of Section VB applies for the quantum harmonic oscillators, as well. 

We now consider two genuine quantum issues, the role of unitarity of the time evolution and of the classical 
limit. First it is pointed out that unitarity protects both reversibility and causality of the dynamics. In fact, both 
irreversibility and acausality are found to be generated by the modified pole structure of the Green-functions, the 
migration of poles to complex location or appearance of new poles. The necessity of poles with non- vanishing imaginary 
part to generate irreversibility indicates that this latter needs non-unitary time evolution. A dissipative force, like the 
friction force in Eq. (87) appears because the large number of environmental soft modes mix with the system state 
and makes it "leaking" to the environment. As of causality is concerned, note first that the final time T of the time 
evolution in the generating functional (71) can be chosen in an arbitrary manner as long as it exceeds the time of any 
observation. In fact, the extension of the time interval the system is followed, T — > T + Ti amounts to the insertion 
of the operator U^(j+;T + T u T)U{-j-;T + Ti,T) within the trace on the right hand side of Eq. (71). But this 
operator is the identity for unitary time evolution, j~ = —j + and does not change the generator functional. In other 
words, any modification of the physical external source j = j + = —j~ is unobservable by looking into expectation 
values at times before the modification, briefly causality can not be lost by a unitary time evolution. 

Decoherence, a necessary conditions for classical limit stands for the suppression of the contribution to the CTP 
path integral for well separated pairs of trajectories. It is advantageous to introduce the Keldysh parameterization 
[8] of the coordinates, x = (x + + x~)/2, x a = x + — x~ and write the O (x 2 ) part of the action in Eq. (89) as 

xD- x x = -xD- la x a + -x a D~ lr x + -x a D- 1 x a , (91) 
2 2 2 y J 

where D~ lr , D~ la and D^ 1 arc defined as in Eq. (56) for crl)" 1 ^, Z) -1 being given by Eqs. (60)-(62). The quadratic 
forms D~ lr {t 1 t') and D~ la (t,t') are real but l)~ 1 (t,t') is imaginary. Hence the magnitude of the contributions to 
the path integral (89) is controlled by the last term. Due to the positive definiteness 3_D" 1 (cj) > 0, l) _1 (i,t') is the 
origin of environment induced decoherence. 

It was pointed out in Section III B that the time arrow of a driven harmonic oscillator is set by the homogeneous 
solution of the equation of motion. In a similar manner we can write the retarded Green-function of the effective 
theory as a sum of a symmetric and asymmetric component, D r = D n + , the latter, being responsible of the 
homogeneous solution and the orientation of the time arrow. Eqs. (66) and (68) show that the time arrow is now set by 
the combination of the spectral functions p + (|u;|) — p_(|w|) and the suppression is controlled by /9 + (|w|)+/9_(|u;|). The 
latter combination leads to decoherence and the former shows the competition between oscillators with oppositely 
rearranged boundary conditions in time in defining the system time arrow. If this combination is non-vanishing 
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then it dominates the infinitesimal D , given by Eqs. (A29) and may generate irreversibility. The lesson is that 
for a usual initial condition problem where po = irreversibility and decoherence stem from the same dynamical 
mechanism. Their difference can be explored mathematically by imposing final rather than initial conditions for some 
environmental degrees of freedom because the irreversibility and decoherence are constructed in the same manner 
except the time arrow of the environment modes are kept or ignored, respectively An extreme manifestation of this 
difference, reversible decoherence is found when the forward or backward moving environment modes have identical 
spectral weight. 

VII. SUMMARY 

The dynamical generation of time arrow of an oscillator is studied in this work within a harmonic model. An 
oscillator has an internal time arrow, defined by its boundary conditions in time, the description of what is supposed 
to be kept fixed when its dynamical environment is changing. The influence of the environment at the normal 
frequency of the oscillator, described by the far Green-function may change the internal time arrow in setting up a 
dynamical time arrow. Irreversibility corresponds to a complete overwriting of the internal time arrow. In cases when 
the time arrow is lost as a result of this process one finds acausality. 

It is proposed that irreversibility and acausality build up in a manner similar to a phase transition. In particular, 
irreversibility may be considered as a dynamical symmetry breaking. This analogy is motivated by the slowing down 
of the order parameter in the thermodynamical limit, the discontinuity of the long time observation limit for infinite 
systems. Such a mathematical subtlety is important if one keeps in mind that observations in physics always deal 
with finite systems and infinity or certain limits seem to be useful approximations rather than part of physical reality 
around us. One may say that the continuous spectrum approximation which allows us to derive irreversibility or 
even the construction of continuum in mathematics by means of the Bolzano- Weierstrass theorem is an ingenious 
approximation which eliminates unreasonable large or small quantities, decoupled from the phenomena of interest. 

The particular mechanism to drive a spontaneous symmetry breaking can easily be modeled phenomenologically 
by introducing an infrared cutoff on the observations. When the dynamics of an infinite system is followed for 
finite amount of time then the extracted laws may reflect less symmetry. In case the frequency spectrum of our 
harmonic model possesses a condensation point with sufficiently strong spectral weight the normal modes around the 
condensation point serve as a sink for the energy received from an external source and observations, carried out during 
finite amount of time signal irreversibility. 

One can always find an order parameter for a spontaneous symmetry breaking. If thermodynamical quantities can 
be defined then the entropy production is a natural order parameter for irreversibility. Entropy is based on limited 
observability, the separation of the full system into an observed subsystem and its environment and only the former 
can be irreversible. This remark opens the possibility that different splitting of the full system into observed and 
ignored parts lead to different status of the time reversal invariance for the observed component. Since dynamical 
breakdown of symmetries requires the cooperation of infinitely many degrees of freedom it is reasonable to expect that 
the ignored component must contain infinitely many sufficiently decorrelated degrees of freedom to show irreversible 
system dynamics. An interesting problem, left open is to find the "irreversible universality class", the subsets of the 
full system which display irreversibility. 

Despite the similarities between irreversibility and spontaneous symmetry breaking these two phenomena show 
important differences, as well. On one hand, the conventional spontaneous symmetry breaking is driven by decorrelated 
domains in three-space. On the other hand, irreversibility happens in time. The dimensionality and the signature 
of the manifold in the relativistically invariant length differ, as well. The dynamical quantities remain correlated 
for arbitrarily large time separation, manifested by conservation laws. Further important non-local effects in time 
are generated in the effective system dynamics by the elimination of degrees of freedom. As a result spontaneous 
symmetry breaking is made possible along the time, in a one-dimensional chain, the less correlated, "disordered" 
phase belonging to irreversible dynamics. A less formal difference is related to initial conditions. Though the common 
dynamical origin of the slowing down of the order parameter and irreversibility is the high frequency level density 
at a system normal frequency these phenomena require different initial conditions. In fact, spontaneous symmetry 
breaking is observed when the energy in the initial state is as low as possible. But irreversibility requires some initial 
energy distributed unequally among the degrees of freedom or the presence of a driving force. The model (1) with 
weak couplings belong to the latter case because the full system starts close to its minimal energy state and energy 
can be gained from the source j(t) only. The energy is received by the system and is dissipated to the environment. 
For sufficiently many soft modes such a dissipation can be maintained for long enough time and the system serves as 
a transit station of energy between the external source and the environment. It is obviously this energy flux which is 
the dynamical origin of irreversibility. There are furthermore aspects, such a locality or local equilibrium which are 
important for phase transitions and can not even be addressed in the model (1). 
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The need of handling infinitely many degrees of freedom to establish irreversibility forces the use of more sophis- 
ticated methods of solving the equations of motion than their direct integration in time which produces trajectories 
with unspecified analytical properties. The algebraic elimination of the environment degrees of freedom in Fourier 
space requires the inversion of some operators which is possible only if their domain, the functional space where they 
act is properly specified. The mathematical background for inversion is easiest to provide within the action principle. 
Thus the CTP formalism is worked out for classical systems where the initial condition problem is transformed into a 
variational problem for periodic trajectories with well defined Green-functions. The near and far Green-functions are 
constructed for a finite observation time and are calculated explicitly in the long time limit where time translation 
invariance can be recovered. 

The reduplication of the degrees of freedom, arising from the periodicity of the trajectories in time is an inherent 
property of the CTP formalism and it explains a surprising feature, the necessity of going beyond the physical space of 
trajectories to establish an action principle for initial conditions which is local in time. This enlargement of the space 
of trajectories has no observable effect in classical mechanics where the equations of motion force the CTP doublers to 
follow the same trajectory. But quantum effects can be encoded in the difference between the doublers owing to the 
linear superposition principle, since the doublers experience independent quantum fluctuations. One may even define 
quantum fluctuations in the path integral formalism by the difference of the doublers trajectories. In the classical 
limit when the system becomes increasingly decohered the doubler trajectories tend to stick together. When the full 
system starts in thermal equilibrium, represented by a canonical density matrix then there is a specially clear way of 
separating thermal and quantum fluctuations as fluctuations in the sum and difference of the doubler trajectories in 
the path integral formalism. 

Decoherence and irreversibility are supposed to be two necessary ingredients of the classical limit. It is found that 
both are determined by the same spectral functions the only difference being that the former is fed by the unoriented 
and the latter by the oriented combination of the spectral functions. 

The advantage of the simple harmonic model considered in this work is its exact solvability. On the other hand, 
obvious difficulties hinder the generalization of the results to more realistic anharmonic models with true interactions. 
Nevertheless, there are good reasons to expect the qualitative results to remain valid for interactive systems, as well. 
Neither irreversibility nor acausality can be generated in any finite order of the perturbation series. Rather, they must 
appear on the quasi-particle level of an interactive many-body system. Our harmonic model correspond to quasi- 
particles obtained by a partial resummation of the perturbation expansion and thus is supposed to be a qualitatively 
acceptable approximation except in the vicinity of the reversible-irreversible transition. It appears desirable to check 
if the generalization to interactive systems, where new time scales arise, confirm our general conclusions. A further, 
more principal question is if the onset of irreversibility is indeed a bona fide phase transition. On one hand, it seems 
to be the simplest realization of the idea of spontaneous symmetry breaking where the initial conditions modify the 
long time dynamics by the slowing down of the order parameter. On the other hand, however spontaneous symmetry 
breaking is a phenomenon where the ground state (respectively the equilibrium state) displays less symmetry than 
the microscopic equations of motion. This is here not the case because the true ground state (or equilibrium state) is 
static by nature and thus does not possess a time arrow. Another, related question is if there exist local parameters 
of the dynamics whose variation drives the system across the reversible-irreversible transition line. This question 
can not even be addressed within our simple model since locality asks for a field theory. A free field model would 
be not sufficient either due to its strongly constrained dynamics. Therefore we are led to the problem of including 
interactions in a field theoretical setting and to search for a nontrivial phase structure with an irreversible phase. 

Another issue raised and left open in this work is the different manner acausality may manifest itself. As usual in 
nonrelativistic physics we define causality as the violation of the order of action and reaction. The question arises 
if the same mechanism would also lead to a spontaneous breakdown of symmetry with respect to Lorentz boosts in 
relativistic models. 
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Appendix A: CTP Green-function 

The formulae for the Green-function of the CTP formalism are collected in this Appendix. First the action is defined 
for fixed IR and UV cutoffs in Section A 1, followed by the derivation of the general expression of the Green- function 
in Section A 2. The explicit formulae of the Green- function for finite and infinite time are presented in Sections A3 
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and A 4, respectively. Section A 5 is for the calculation of the inverse Green-function, the quadratic form of the kinetic 
energy in the infinite time case. 

1. Action 

The identification of the end point of the trajectories x+ and x~ requires care and we consider a harmonic oscillator 
with mass m and frequency w for discrete time, t = jAt, j = 0, 1, . . . , Nuv, At = T jNjjy. The real part of the CTP 
action is 



MS, 



CTP 



m 
2At 



t+At 



x lf 



AtnvJk 



E 



X, 



(Al) 



t t 

where the sum over t is for < t < (Nuv — I)At. Partial integration/summation in the kinetic energy yields the form 



StScTP — E " Ta"7 E X t(^ X t X t-At X t+At) + ^~L[ x NAt( x NAt — x (N-l)At) x o( x At 



X 



AtmujQ 
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The constraint 



S c 



1 1 v , I . 

^Pi — At At XAt ' 



handling the initial condition gives the extended action 

sctp - 1 E E *t D o\y'*t +EE^( A " + 



o,tj' t,t> 

where the notation y = x^At has been introduced and 

-lCT.cr' 



<T t 



D 



d t.V 



a ( A7^*'*' + ^ tmuj o s t,t>^J - iAtmedt. 



Ira 
At ' 



(A3) 



(A4) 



(A5) 



This expression contains an UV and an IR cutoff, At and T. The UV cutoff is removed first in the continuum limit, 
At — > by keeping T fixed. The way the IR cutoff is treated is determined by our choice of the scheme NT or TN. 
The removal of the cutoffs is nontrivial because neither of the limits At — > or T — >• oo is continuous. 



Green-function 



To find the Green- function for the quadratic action (53), (A4) we write 

Sctp = \xD^x + x(A + Bz + j). 



(A6) 



We keep At finite which allows to make the trajectories x + (t) and x~(t), given for < t < T independent by treating 
the final point, z — x ± (T) separately. Due to the absence of O (z 2 ) term in the action we eliminate first the trajectory 
x by means of its equation of motion, 



x = -D d (A + Bz+j) 
and express the action in terms of the variable z only, 

Sctp = ~(A + zB + j)D d (A + Bz+j). 



(A7) 



(A8) 



24 



The variational equation for z, 



z = _BDM + j) 
BD d B 



gives finally the action as a functional of the source, 



by means of the free Green-function 



W\j] = ~(A + j)D (A + j), (A10) 



b Q = D d -b d B^-L^-BD d . (All) 
BD d B 



The initial condition can be imposed on the trajectory 



x = ^- = -D (A + j) (A12) 



by the appropriate choice of the multiplicative factor t in A. 

3. Finite T 

To find the inverse of I) 1 for Eq. (All) we have to specify the functional space where we are working. It is span 
by functions x(t) defined over the time interval [0, T] and satisfying the boundary conditions x(0) = x(T) = 0, 

Nuv IY 

x(t) = ^2x n J— smu n t, (A13) 

n=l 

with u) n — nir/T and non-uniform convergence is allowed at t = T in the continuum limit, Nuv — > 00 to reach an 
arbitrary end point z = lim t ^.T x(t). 

The starting point is the Green-function 

Dd (t, t = — > — 2—— A14 

where u) n = sin7rA|2.. The continuum limit has to be taken with some care. The usual procedure is to split the 
right hand side into the sum for 1 < n < cNjjv and for cNjjv < n < Nuv where c is a small, TV^y-independcnt 
number. The latter sum, being proportional to l/TVyyc 2 can be neglected and the replacement £b n — > uj n can be made 
in the denominator and we get the diagonal CTP block-matrix 

D d - ft _" D .) (A15) 



with 



n u +"\ 2 ^ sin uj n t sin uj n t' 

D ^^ = ^2^ (A16) 



after the removal of the UV cutoff. 

We need in the denominator in the right hand side of Eq. (All) 

1 D(T At T At) 2 y shi^Cr -At) 

n— 1 " u 

that can be split into the sum of an O {Nuv) linearly divergent part and a finite piece, 



-^D d (T — At, T - At) = —— lim V 1 - sin 2 ^ + ±% V 9 2 . , (A18) 
At 2 mTN uv -*oo^\ 2 J mT ^ w 2 - uj 2 + it y ' 
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in the continuum limit. Another ingredient of Eq. (All) becomes in the limit Nuv —> oo 



1 „ , „ a \ 2 lo,„ sin ui n (t — T) , . , 



At ° v ' ' Tm^ w^-w^+ie 

n=l n u 



The complete Green- function (All), 

ft ( . fD d (t,t>) \ im ( D d {t,T - At)D d {T - At,t<) D d (t, T - At) D* (T - At, f)\ 

0( ' >-{ -D*(t, t') ) 2 zw 2 At 2 ^3(*. T At)D+(T - At, t') D* d (t, T - At)^(T - At, t')J [AZ[)) 



with 



converges in the continuum limit. 



4. Diverging T 

The Green- function can directly be determined for T = oo as argued in Section VA where it was noted that the 
limit T — > oo is discontinuous. Our goal here is to find the modification of the rules of the limit T — > oo to recover 
the result for T = oo. 

First we execute the limit T — > oo for fixed small but finite e. The diagonal Green-function, D d , given by (A16) is 

DM, t') = !_(e-* 5 °l*-*'l - e -^o(t+0). ( A22 ) 

2mwo 

with wo = Wo — ie/2u>o- The expression (A18) becomes 

1 n 1 .wo 

— o ^OT-At.T-At — 7— - « 

At mAt m 

in the limit At — > 0. The expression (A19) requires more care in order to correspond to the motion between the time 
t and T. One finds 

^-D (t, T-dt) = — L e -«>o( r -*) (A24) 
At m 

by keeping t — T < fixed in the limit T — > oo. We expect translation invariance to recover far from both the initial 
and the final time hence we seek the Green-function by making the replacement t — > T/2 + t with fixed t as T — > oo. 
The substitution of these results into eq. (All) yields 

'g— iwol*— 1'| o 



2 D 0T -At,T-At = ^7 - i— (A23) 
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e iui (t+t' -T) e iui (t-t') \ / g— iw (t+t') 
g-iwo(t-t') e -iw (t+t'-T) J _ I q e iw (t+t') 



)]}' 



(A25) 



A simple form, given by Eq. (59) reproduces Eqs. (28) and yields 



- , , cos Wot 

A>(* = L - A2 6 

This Green-function (59) can be recovered if the limit T — > oo is performed by first approximating the infinite sums 
in Eqs. (A16), (A18), (A19) and (A21) for large but finite T by integrals and after that by performing the limits 
e — > and T — > oo in such joint manner that the limits 

e "2^ T i (A27) 

apply within the integrand of the right hand side of (11). The two steps contradict because eT — > is required by 
the second limit in (A27). However the infinite sum can not be approximated by a contour integral if the distance 
of the pole from the real axes is smaller than the discrete increment of the summation variable. In other words, the 
continuous spectrum formalism can not be obtained as the infinitely long time limit. The formal solution is based on 
contour integrals where an infinitesimal e is introduced which is strong enough for the first limit in (A27) and weak 
enough for the second. This is in agreement with the practice of using ie as a prescription to avoid poles rather than 
a finite quantity. 
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5. Inverse Green-function 



The translation invariant Green-function (59) can be used to find the CTP action for T = oo. Its Fourier transform, 

D ^ = „ l-Mewfa - ug) J ' 

shows the structure mentioned in Section V A, namely the trajectory (54) is constructed by collecting the unoriented 
driven modes by D n (ui) from one time axis and the free modes are taken from the other time axis by D*(u>) in such 
a manner that the sum is a retarded solution. 

To find the inverse of Dq{u)) we use S(x) = e/ir(x 2 + e 2 ) for the regulated Dirac-delta. The trivial algebraic inversion 
yields the form (60) or 

D Q ln = m(u? - uj 2 ), 
Dq 1 ^ — isign(fc°)e, 

£>o" = e - ( A29 ) 

The coupling of the two trajectories which is at the final point on the original action (50) is hidden in the limit T — > oo 
and is replaced by an infinitesimally strong coupling distributed along the whole time axis in a uniform, translation 
invariant manner. 
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